Effect of electron irradiation on vortex dynamics in YBa2Cu307_5 single crystals 
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We report on drastic change of vortex dynamics with increase of quenched disorder: for rather weak 
disorder we found a single vortex creep regime, which we attribute to a Bragg-glass phase, while for 
enhanced disorder we found an increase of both the depinning current and activation energy with 
magnetic field, which we attribute to entangled vortex phase. We also found that introduction of 
additional defects always increases the depinning current, but it increases activation energy only for 
elastic vortex creep, while it decreases activation energy for plastic vortex creep. 



PACS numbers; 74.60. Jg, 74.60. Ge, 74.60. Ec, 74.72.Bk 



The effect of random pinning on crystalline order and 
on dynamics of the flux-hne-lattice (FLL) was a subject 
of numerous experimental and theoretical investigations. 
Neutron diffraction § and /iSR §1 experiments gave ex- 
perimental evidence for the existence of two vortex solid 
phases with different positional correlations. Magnetiza- 
tion measurements ||] of Bi2Sr2CaCu20s showed that 
the field, Hon, at which a steep increase of the measured 
current Jm{H) begins, coincides approximately with the 
field above which the intensity of Bragg peaks sharply 
decreases The field Hon was interpreted as a phase 
boundary between low field ordered and high field dis- 
ordered vortex phases. These experimental results are 
supported by theoretical studies. It was shown that in 
presence of rather weak disorder the FLL retains a quasi- 
long-range order resulting in the so-called Bragg glass 
phase However, with the increase of random pin- 

ning or magnetic field a transition to strongly disordered 
entangled vortex phase (glass phase) is predicted . A 
sharp increase in magnetization below the fish-tail peak 
position Hp was observed in Ndi,sbCeo,ibCu04^-d @] 
and non twinned YBa2Cu307_5 [D single crystals. 

Magnetic measurements showed that in optimally 
doped YBa2Cu307_5 crystals no fishtail behavior is ob- 
served both in detwinned Q and twinned samples, 
while decrease of the oxygen content always induced non 
monotonous Jm(-ff)-curves. Two distinctive peculiarities 
in the magnetization curves of oxygen deficient | |To[ | and 
electron irradiated |11| crystals were observed: (1) the 
peaks Hon and Hp shift toward lower magnetic fields with 
increasing defect concentration, and (2) in magnetic fields 
increases, while in magnetic fields 
decreases with increasing defect 
concentration. It is believed [ p^|p^ that the peak Hp 
separates elastic vortex creep in low and plastic vortex 
creep, mediated by motion of the FLL dislocations, in 



H < Hp the current J„ 
H > Hp the current J, 



high magnetic fields. Thus the introduction of additional 
defects leads to an increase of vortex pinning in the re- 
gion of elastic creep and such behavior is expected. On 
the other hand, decrease of pinning force with increas- 
ing disorder observed in the region of plastic creep is non 
trivial, and reasons of such behavior are not known. The 
aim of this paper is to show the effect of point-like defects 
concentration on vortex dynamics and pinning parame- 
ters in YBa2Cu307_5 single crystals. 

The investigated sample was YBa2Cu307_5 single 
crystal with « 93.5 K and AT^ < 0.5 K. Twin 
planes (TP's) inside the measured part of the sample 
were aligned in one direction. The transport current was 
applied along the a5-plane and at an angle a = 45° with 
respect to TP's. Measurements were performed in mag- 
netic fields applied parallel to the c-axis. Temperature 
stability of the sample during measurements was bet- 
ter than 5 mK, and measurements in the normal state 
showed that overheating of the sample at the highest dis- 
sipation level of 50 /iW did not exceed 10 mK. 

Additional defects were introduced by irradiation with 
2.5-MeV electrons, which are suitable for production of 
point-like defects. Irradiation was performed at temper- 
atures T < lOK and after irradiation the current- 
voltage characteristics (CVC) were measured without 
heating the sample above 110 K. This excludes diffu- 
sion, and therefore annihilation and clustering of the 
defects. Irradiation was performed at dose rate 4.2 x 
10^^cm~^sec~^ and at an angle 5° off the c-axis to avoid 
electron channeling. Homogeneity of electron beam was 
about 5%. Following procedure used in Ref. (l5| we esti- 



mated that irradiation dose IO^^cto"^ produces the aver- 
aged over all sublattices concentration of the defects 10~^ 
dpa, and that the penetration range ~ 1mm is at least 
two orders higher then thickness of the crystal = 7/xm 
resulting in homogeneity of the defects along the c-axis 
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better then 1%. Reduction of the Tn after irradiation was 



0.65, 1 and 1.3 K after irradiation with doses 



10^^,2 X 10^^ and 3 x lO^^e^/cm^, respectively. 

Fig.l shows field variation of the measured "criti- 
cal" current determined at an electric field level 
10~^V/cm. Before irradiation the current J,„ continu- 
ously decreases with increasing field as it was previously 
observed by Zhukov et al. in crystals with very small 
oxygen deficiency, 5 ~ 0.03 ||l^]. After irradiation the 
Jm{H) curves show fish-tail behavior and the peak po- 
sition Hp in the Jm{H) curves shifts toward low fields 
with the dose. It is also seen that introduction of addi- 
tional defects increases Jm for H < Hp , while the current 
Jm decreases with the dose for H > Hp. Such variation 
of the Jm and Hp with increasing of the defect concen- 
tration is analogous to the behavior previously observed 
in oxygen deficient |]l0| and in electron irradiated ||l|] 
YBa2Cu307_5 single crystals. 
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FIG. 1. Field variation of the current Jm- The inset shows 
vortex velocity versus J before (open symbols) and after (dark 
symbols) irradiation with (p — 10^* el/cm^. 

In low magnetic fields {H < 15kOe for non irradiated 
sample, and H < Hp for irradiated samples) the CVC 
data follow the dependence 



E^Eoexp[-{Uo/kT){JjJY 



(1) 



where the exponent = 1, i?o is a constant, Uq is 
the activation energy, and k is the Boltzmann constant. 
The depinning critical current Jc can be determined by 
extrapolation of the ratio Pd{J) / Pbs to unity fl^ as- 
suming that at current density J = Jc the differen- 
tial resistivity pd = dE/dJ equals the flux flow resis- 
tivity in the Bardeen-Stephen model, pss — PnB/Bc2 
where pat is the normal state resistivity. Induction 
of the upper critical field was estimated assuming that 
Bc2 = {dBc2/dT){T -Tc) wfth dBc2/dT = -2.5T/K M. 
Field variation of the Jc is shown in the inset of Fig]|a. 
Substituting these values of Jc in Eqj|, and fitting exper- 



imental £'(J)-curves by this equation we derived Uq{H) 
dependence shown in the inset of Fig.^ja. 
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FIG. 2. 
The inset 
The E{J) 

radiation with doses 



(a) The E{J) curves plotted as E{J)vs.l/J. 
shows field variation of the Jc and Uo. (b) 
curves measured before (squares) and after ir- 
(circles), 2 x 10^** (up triangu- 



10^ 



lars), and 3 x 10 el /cm (down triangulars) , and plotted as 
E{J) / Jvs.J^^^ . The inset shows Jp; (_ff)-dependence (dark 
symbols) and (/p; (-ff )-dependence (open symbols) derived for 
T = 85K before irradiation (squares) and after irradiation 
with doses 10^* (circles), 2 x 10^* (up triangulars), and 
3 X 10^^ el /cm? (down triangulars), and derived for T = W.K 
and = 3 X IQ^^el/cm^ (diamonds). 



As one can see in Fig.|l| and in the inset of Fig.|^, vor- 
tex dynamics, namely field variation of the pinning pa- 
rameters, strongly depends on the strength of disorder. 
In presence of weak disorder (before irradiation) Jc^ Uo, 
and vortex velocity v — cE / B do not depend on mag- 
netic field indicating a single vortex creep regime, as it 
is predicted by the collective pinning theory |16| for low 
magnetic fields. Also, for this regime of vortex creep the 
current J„ decreases with increasing magnetic field due 
to increase of the fiux density. Increase of the disorder 
induces field variation of the pinning parameters. After 
irradiation Jc and Uq increase with magnetic field that 
leads to reduction of vortex velocity and to increase of 
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the current Jm with increasing magnetic field. These ob- 
servations indicate that before and after irradiation we 
test different vortex phases. 

In presence of weak quenched disorder the ordered 
Bragg-glass phase is expected Q and our experimental 
data probably indicates that dynamics of just this vortex 
phase is described by the single vortex creep regime. Also 
note, that derived exponent /x = l±0.15is close to the 
value of /i = 0.7 — 0.8 predicted for the Bragg-glass phase 
for not too small currents Q . With increase of disorder 
a transition of the Bragg-glass to entangled vortex phase 
is expected |^,^- The derived after irradiation weak in- 
crease of the Uq and rapid increase of the Jm with mag- 
netic field correlates with previous experimental findings 
for magnetic fields B > Bon, or in magnetic fields 
where the entangled vortex solid is expected. Therefore 
we believe that after irradiation we test the entangled 
vortex phase and dynamics of this phase is characterized 
by the increase of both Jc and Uq with increasing mag- 
netic field, which lead to rapid increase of the current 
J„i. Also, introduction of additional defects increases 
both the Uq and Jc resulting in steep increase of the cur- 
rent Jm with irradiation dose, see Fig.|l|, in agreement 



with previous investigations |11 



Let us consider experimental data obtained in mag- 
netic fields H > 20kOe for nonirradiated sample, and in 
magnetic fields H > Hp for irradeated samples. As one 
can see in Fig.^, the CVC data follow the dependence 
predicted for motion of the FLL dislocations MJi^ 



EiJ) = poJexp{-iUpi/kT)[l - [Jpi/JT]}, 



(2) 



where = -1/2, and po is a constant. The main pecu- 
liarities of the curves presented in Figj^b are that the 
slope of the curves measured in the same magnetic field 
decreases with increasing defects concentration and that 
these curves intersect one another. For vortex creep de- 
scribed by Eq.^ such behavior is possible only in that case 
when the critical current increases, but the activation en- 
ergy decreases with increasing defect concentration. Ex- 
trapolating the ratio Pd{J) / Pbs to unity we derived field 
and temperature variation of_the current Jpi shown in 



2 a, respectively. Substitut- 



the inset of Fig.pD and in Fig, 
ing these values of Jpi in Eq. 2 , and fitting experimental 
i?( J)-curves by this equation we determined the field and 
temperature variation of Upi shown in the inset of Fig.||b 
and in Fig.^, respectively. 

The current Jpi is determined by the interaction of 
the dislocations with the FLL and random pinning cen- 
ters. Interaction with the FLL results in a shear lim- 
ited contribution |^ Jgh « cee/Bd oc i?-^/^, where 
C66 = <i>o-B/(87rA)^ is the shear modulus, $o is the flux 
quantum, A is the penetration depth, d « oq is the width 
of channel for moving dislocation, and gq « (^o/B)^^^ 
is the intervortex distance. As one can see in the inset 
of Fig.||b, in magnetic fields not very close to the melt- 
ing point the current Jpi really increases with the field 



in agreement with theoretical predictions. In the pres- 
ence of random pinning the shear modulus cge decreases 
PH and therefore Jsh also decreases. However, the con- 
tribution of the core pinning increases with the defect 
concentration . The derived increase of the current 
Jpi with irradiation dose indicates that increase of the 
core pinning dominates over reduction of the Jsh- 
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FIG. 3. (a) Temperature variation of Upi and Jpi for 
H = 15fcOe. (b) Temperature variation of Jm for H — IbkOe. 
The inset shows nuclei for the plastic creep in the absence 
(left) and in the presence (right) of random pinning. 

The activation energy Upi decreases with the increase 
of both irradiation dose and magnetic field. A decrease 
of the activation energy with increasing field agrees with 
theoretical calculations, Upi « eeooo oc H, where 

£o = (*I*o/87rA)^, and previous experimental findings 
[|I2|,^. But, a decrease of the activation energy with 
increasing irradiation dose is not evident. Of course, the 
Tc decreases after irradiation which leads to an increase 
of the A(T) = A(0)/[1 - (T/Tc)^]^/^, and hence to a de- 
crease of the activation energy Upi cx A~^. However, the 
reduction of the Tc can not describe the fast decrease of 
Upi in our measurements. This is demonstrated in FigJ|a, 
which shows dependence Upi vs. t = T/ (Tc — ATirr)- It is 
evident that Uq decreases with the dose increasing even 
if we take into account the reduction of T^. 

To explain this behavior let us consider a displacement 
of a vortex segment Lq over intervortex distance oq as 
shown in the left-hand inset of FigJ|b. In the absence of 
random pinning the energy of such nucleus can be writ- 
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ten as Upi = 2Eei + Epi, where Eei ~ eeooo is elastic 
energy of the vortex segment Lab^ and Epi is the energy 
required for displacement of the vortex segment over 
the distance oq. At small driving forces such a nucleus 
is stable when Epi > 2Eei. Thus we obtain a minimal 
activation energy Upi ~ -ieeoao, which coincides within a 
factor 4 with previous estimates Upi w eeooo In the 
absence of random pinning the equilibrium positions of 
the vortices are the straight lines as it is shown by two 
parallel straight lines. In the presence of random pin- 
ning the equilibrium positions become curved as shown 
by dashed lines in the right-hand inset of Fig.||b. There- 
fore, along the vortex lines a vortex fragments, for which 
an average distance Up between two neighboring equi- 
librium positions of vortex lines is smaller then qq ap- 
pear. In this case the energy Upi w iesoUp is smaller 
then activation energy in the absence of pinning because 
flp < flQ. Substantial displacements of vortex segments, 
Aao = oo — flp ~ 0.2ao ~ 8nm, due to core interaction 
with individual pinning centers probably unreliable, and 
we attribute them to fluctuations of the defects concen- 
tration, which naturally present in real crystals. Indeed, 
equating the work A = JpiFoLoUp/c, required for dis- 
placement of the segment Lq over distance Op, to the elas- 
tic energy, 2Eei = Upi/2, we estimated Lq = 500 nm for 
ip = 10^^ el/cm^, and derived that for coherence length 
^(85K) = 5 nm the segment Lq interacts with about 200 
point defects. Therefore, fluctuations of the defects con- 
centration of 10% can give difference in core interaction 
with about 20 point defects in different positions of the 
segment Lq. 

TPs strongly affect pinning and dynamics of vortices 
p3|"p5[. In particular, being plane defects they form 
channels of easy vortex motion along the plane of twins 
P^ , p5| . Due to suppression of superconducting order pa- 
rameter within the TPs some part of vortices is trapped 
by the TPs [Q, and pinning of these vortices along the 
TPs may be reduced as compared with pinning in the 
bulk of the crystal ||2^. Therefore, for the same driving 
force velocity of the trapped vortices can be higher com- 
pared with velocity of vortices placed in the bulk of crys- 
tal. In high magnetic fields contribution of the trapped 
vortices to dissipation of energy is small due to small frac- 
tion of these vortices. However, in a magnetic field 1 kOe 
the intervortex separation uq = 140nm becomes compa- 
rable with the distance between twins d = SOOnm in our 
sample, and a significant part of vortices can be trapped 
by the TP's. Therefore, contribution of the trapped vor- 
tices increases with decreasing magnetic field, and devi- 
ation from the field scaling for non irradiated sample in 
low fields presented in the inset of Fig jl] probably reflects 
the reduced pinning of the trapped vortices. 

In conclusion, we have shown that vortex dynamics 
dramatically depends on the strength of disorder. In 
presence of weak disorder a single vortex creep is real- 
ized, which we attribute to dynamical property of the 



ordered Bragg-glass phase. In presence of strong disorder 
and for elastic vortex creep we have found rapid increase 
of the depinning current and weak increase of the acti- 
vation energy with increasing magnetic field, which we 
attribute to dynamical property of the entangled vortex 
solid. We also found that the introduction of additional 
defects always increases the depinning critical current, 
but it increases the activation energy only for the elastic 
vortex creep, while it decreases the activation energy for 
the plastic vortex creep. 
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